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Abstract 

We calculate the one-photon loop radiative corrections to virtual pion Compton scat- 
tering 7r~7* — > vr~7, that subprocess which determines in the one-photon exchange ap- 
00 ■ proximation the pion-nucleus bremsstrahlung reaction ■k~ Z ^ 'k~ Z^. Ultraviolet and 

\ infrared divergencies of the loop integrals are both treated by dimensional regularization. 

■ Analytical expressions for the 0{a) corrections to the virtual Compton scattering am- 
plitudes, ^(s,M,Q) and B{s,u,Q), are derived with their full dependence on the (small) 

Q I photon virtuality Q from 9 classes of contributing one- loop diagrams. Infrared finiteness of 

' these virtual radiative corrections is achieved (in the standard way) by including soft pho- 

f-*) ■ ton radiation below an energy cut-off A. In the region of low 7r~7 center-of-mass energies, 

where the pion-nucleus bremsstrahlung process is used to extract the pion electric and 
magnetic polarizabilities, we find radiative corrections up to about —3% for A = 5 MeV. 

■ Furthermore, we extend our calculation of the radiative corrections to virtual pion Comp- 
ton scattering 7r~7* tt^j by including the leading pion structure effect in form of the 

Oh| polarizability difference a-,^ — [i-K- Our analytical results are particularly relevant for ana- 

lyzing the data of the COMPASS experiment at CERN which aims at measuring the pion 
electric and magnetic polarizabilities with high statistics using the Primakoff effect. 



> 

cn . PACS: 12.20.-m, 12.20.Ds, 13.40.Ks, 14.70.Bh 

1 Introduction and summary 

00 ■ 

^ ' At present there is much interest in a precise experimental determination of the pion electric 
' and magnetic polarizabilities, a^r and (3-,^. Within the framework of current algebra [1] it has 
^ ■ been shown (long ago) that the polarizability difference a-,, — P-,, of the charged pion is directly 
^ . related to the axial-vector-to-vector form factor ratio hA/hy ~ 0.44 measured in the radiative 
pion decay 7r+ e^Vel [2]- At leading (nontrivial) order the result of chiral perturbation 
theory [3], — = a{lQ — r5)/(247r^/^m^), is of course the same after identifying the low- 
energy constant as Iq — h = GhA/hy- Recently, the systematic corrections to this current 
algebra result have been worked out in refs.[4, 5] by performing a full two-loop calculation of 
(real) pion Compton scattering in chiral perturbation theory. The outcome of that extensive 
analysis is that altogether the higher order corrections are rather small and the value — /^jr = 
(5.7 ± 1.0) ■ 10~^fm^ [5] for the pion polarizability difference stands now as a firm prediction of 
the (chiral-invariant) theory. The non- vanishing value a,r + /^tt = (0.16 ±0.1) ■ 10~^fm^ [5] for 
the pion polarizability sum (obtained also at two-loop order) is presumably much too small to 
cause an observable effect in low-energy pion Compton scattering. 

However, the chiral prediction — /^^r = (5.7 ± 1.0) ■ 10~^fm^ is in conflict with the existing 
experimental determinations of — = (15.6 ± 7.8) ■ 10~^fm'^ from Serpukhov [6] and 
— = (11.6 ± 3.4) ■ 10~'^fm^ from Mainz [7], which amount to values more than twice as 
large. Certainly, these existing experimental determinations of — /Stt raise doubts about their 
correctness since they violate the chiral low-energy theorem notably by a factor 2. 
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In that contradictory situation it is promising that the ongoing COMPASS [8] experiment 
at CERN aims at measuring the pion polarizabihties with high statistics using the Primakoff 
effect. The scattering of high-energy negative pions in the Coulomb field of a heavy nucleus 
(of charge Z) gives access to cross sections for 7r~7 reactions through the equivalent photons 
method [9]. In practice, one preferentially analyzes the spectrum of bremsstrahlung photons 
produced in the reaction -k^ Z ^ -k^ in the so-called Coulomb peak. This kincmatical regime 
is characterized by very small momentum transfers to the nuclear target such that virtual pion 
Compton scattering 7r~7* — > 7r~7 occurs as the dominant subprocess (in the one- photon ex- 
change approximation). The deviations of the measured spectra (at low 7r~7 center-of-mass 
energies) from those of a point-like pion are then attributed to the pion structure as repre- 
sented by its electric and magnetic polarizability [6] (taking often + /^tt = as a constraint). 
It should be stressed here that the systematic treatment of virtual pion Compton scattering 
7r~7* — > 7r~7 in chiral perturbation theory yields at the same order as the polarizability differ- 
ence ctTT — /^TT a further pion-structure effect in form of a unique pion- loop correction [10, 11] 
(interpretable as photon scattering off the "pion-cloud around the pion"). In the case of real 
pion Compton scattering 7r~7 — > 7r~7 this loop correction compensates partly the effects from 
the pion polarizability difference [12]. A minimal requirement for improving future analyses of 
pion-nucleus bremsstrahlung 7r~Z — > Z^ is therefore to include the loop correction predicted 
by chiral perturbation theory. The relevant analytical formulas for doing this have been written 
down in the appendix of ref.[ll]. 

At the required level of accuracy it is also necessary to include higher order electromagnetic 
corrections to the pion-nucleus bremsstrahlung reaction 7r~Z — > 7r~Z7. Such radiative correc- 
tions have been considered some time ago in the works of Akhundov et al. [13] and they have 
been implemented into the numerical program called RCFORGV. Unfortunately, no accessi- 
ble sources to the underlying analytical expressions for the photon-loop amplitudes (which are 
necessary for an independent implementation into new data analyses) have been given by these 
authors. The purpose of the present work is to fill this gap. We will re-evaluate the one-photon 
loop radiative corrections to virtual pion Compton scattering 7r~7* — > 7r~7 and present for each 
class of diagrams the corresponding analytical expressions. The calculation gets considerably 
simplified by considering the T-matrix of 7r~7* — » 7r^7 in the laboratory frame and exploiting 
the circumstance that the time-like polarization vector (1,0) of the virtual photon is orthogonal 
to the space-like polarization vector (0, e*) of the real bremsstrahlung photon. Incorporating 
furthermore crossing-symmetry, one finds that only two dimensionless amplitudes, A[s,u,Q) 
and S(s, Q), depending on three dimensionless Lorentz-invariant kinematical variables, need 
to be specified. 

Our paper is organized as follows. In section 2, we review the cross section for pion-nucleus 
bremsstrahlung (in the one-photon exchange approximation) and define the T-matrix of virtual 
pion Compton scattering tt~j* — > 7r~7 in the laboratory frame together with the tree-level 
Born contributions. Section 3 is devoted to the calculation of the one-photon loop radiative 
corrections. We present analytical expressions for the amplitudes A{s,u,Q) and B{s,u,Q) of 
order a as they emerge from 9 classes of one-photon loop diagrams. Dimensional regularization 
is used to treat both ultraviolet and infrared divergencies of the loop integrals. While the 
ultraviolet divergencies drop out in the renormalizable scalar quantum electrodynamics at work, 
the infrared divergencies get removed (in the standard way) by including soft photon radiation 
below an energy cut-off A. We evaluate in section 4 the resulting finite (real) radiative correction 
factor with a small energy cut-off A introduced in the 7r~7 center-of-mass frame. Numerical 
results are presented for this part of the radiative corrections in the region of low 7r~7 invariant 
masses, where pion-nucleus bremsstrahlung is used to extract the pion polarizabihties. The 
detailed evaluation of the finite (virtual) radiative correction factor is not possible at this stage. 
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since it requires specification and implementation of the actual experimental conditions (such as 
pion beam energy, constraints on the detectable energy and angular ranges, etc.). In sections, 
we go beyond the works of Akhundov et al. [13] and extend the calculation of the radiative 
corrections to virtual pion Compton scattering 7r~7* 7r~7 by including the leading pion 
structure effect (at low energies) in form of the polarizability difference a,^ — In an effective 
field theory approach this feature can be elegantly handled by a suitable (gauge-invariant) two- 
photon contact- vertex. In the appendix, we study some two-photon exchange processes specific 
to scalar quantum electrodynamics which turn out to be negligibly small. 

Our analytical results for the radiative corrections to virtual pion Compton scattering can be 
utilized for analyzing the data of the COMPASS experiment at CERN which aims at measuring 
the pion polarizabilities via high-energy pion- nucleus bremsstrahlung t:~Z ^ t:~Zj. 



2 Cross section for pion-nucleus bremsstrahlung 

In order to keep this paper self-contained, we start with reviewing the cross section for pion- 
nucleus bremsstrahlung in the one-photon exchange approximation. Consider the process, 
'^"{Pi) + ^rcst — ^ 7r~(p2) + lik^e) + ^recoil, whcrc a virtual photon transfers the (small) mo- 
mentum q = p2 + k — pi from the nucleus (of charge Z) at rest to the (high-energetic) charged 
pion. The extremely small recoil energy —go = (f^/SMnuci ~ of the nucleus can be (safely) 
neglected. The corresponding fivefold differential cross section in the laboratory frame reads 
[14]: 



H (1) 
diodQ,rydQ,T^ 'K'^\pi\\q\'^ ' 

with a = 1/137.036 the fine-structure constant and H a squared amplitude summed over the 

— * 

final state (real) photon polarizations {e - k — 0): 



H = 



\Api sin 6^1 p + |5p2sin6'2p + 2Re{AB*)\pi\\p2\ sin sin 6*2 cos 0. (2) 

It has been derived from a T-matrix for the virtual Compton scattering subprocess, 7r~7o — >■ 
7r~7, of the form Tiab = —S7Ta{Ae* ■ pi + Be* ■ P2) with A and B two (in general complex) 
amplitudes. Since the electromagnetic interaction with the heavy nucleus involves only the 
charge density, the polarization four- vector of the virtual (Coulomb) photon has only a time- 
component: (1,0). The quantities 9i and 62 in eq.(2) are the angles between the momentum 
vectors k and pi and between k and p2) respectively. Moreover, (/> is the angle between the two 
planes spanned by {k,pi) and {k,p2). The photon energy is denoted hj uj = \k \ = Ei — E2, 

where Ei = \J\pi\'^ + ml and E2 = y^|p2p + m^ are the initial and final pion energies in the 
laboratory frame. 

Taking into account crossing-symmetry, the T-matrix for virtual pion Compton scattering, 
'^~{Pi) + 7o(?) ~^ '^'{Pi) + e), in the laboratory frame can be further decomposed as: 

Tiab = -Pi \e^ A{s, u, Q) + E2 B{s, u, Q)] + e* \e^ B{u, s, Q) + E2 A{u, s, Q)] \ . (3) 

Here, we have introduced the (convenient) dimensionless Mandclstam variables: s = {p2 + 
ky/ml, u = {pi — ky/m1, t = (pi —p2Y/ml and Q = —q^/m^ (given by Lorentz-squares 
of momentum four- vectors) , which obey the constraint s-\-t-\-u — 2 — Q. In the physical 
region one has: s>l,M<l,t<0 and Q > 0. The less obvious condition u — 1 < Q 
can be deduced by considering its left hand side in the 7r~7 center-of-mass frame. Due to the 
prefactor STra/m^ in eq.(3) the two independent amplitudes A{s,u,Q) and B{s,u,Q) become 
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Figure 1: Tree diagrams of scalar QED for virtual pion Compton scattering. The virtual photon 
on the left hand side couples to heavy nucleus. 

also dimensionless quantities. The crossing symmetry of the T-matrix, which corresponds to an 
invariance under the transformation {Ei,pi) {—E2, —P2) and s m, is manifest according 
to the decomposition written in eq.(3). 

The advantage of working with dimensionless amplitudes and variables shows up already 
when calculating the Born terms for a point-like pion. The (three) tree diagrams of scalar 
quantum electrodynamics shown in Fig. 1 lead to the following extremely simple expressions: 



A{s,u) 



(tree) 



0. 



B{s,u) 



(tree) 



U 



(4) 



where we suppress from now on the dependence of A{s, u, Q) and B{s, u, Q) on the third variable 
Q. Note that the right diagram in Fig. 1 involving the two-photon contact- vertex of scalar QED 
vanishes, due to the orthogonality of the time-like virtual photon and the space-like real photon 
polarizations, eo = 0. The s-channel pole term (generated by the left diagram in Fig. 1) enters 
into the T-matrix through the (crossed) amplitude !?(«, s)^*^*^*^-* = 2/(s — 1). The contributions 
of the pion polarizabilities (predominantly their difference — /J^r — —2/3^) to the amplitudes 
A{s,u) and B{s,u) can be found in eq.(27). 



3 Evaluation of one-photon loop diagrams 

In this section, we present analytical results for the one-photon loop radiative corrections of 
order a to the invariant amplitudes A[s, u) and B{s, u) of virtual pion Compton scattering. We 
use the method of dimensional regularization to treat both ultraviolet and infrared divergencies 
(where the latter are caused by the masslessness of the photon). The method consists in 
calculating loop integrals in d spacetime dimensions and expanding the results around d = A. 
Divergent pieces of one- loop integrals generically show up in form of the composite constant: 

containing a simple pole at = 4. In addition, 7^; = 0.5772 ... is the Euler-Mascheroni number 
and /i an arbitrary mass scale introduced in dimensional regularization in order to keep the 
mass dimension of the loop integrals independent of d. Ultraviolet (UV) and infrared (IR) 
divergencies are distinguished by the feature of whether the condition for convergence of the 
(i-dimensional integral is d < 4 or c? > 4. We discriminate them in the notation by putting 
appropriate subscripts, i.e. ^uy and ^/r. The basic aspects of mass (m^) and coupling constant 
(e) renormalization in scalar QED have been discussed in section 2 of ref. [11]. 

The gauge condition eo = eliminates those (five) loop diagrams where both the incoming 
virtual photon and the outgoing real photon are connected by the two-photon contact-vertex 
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Class: I 
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:> 



Figure 2: One-photon loop diagrams of class I. The virtual photon comes in from the left. 



Class: II 



:> 



L 

^^^^ 



Figure 3: One-photon loop diagrams of class II. 



of scalar QED (in first or second order). One can also drop the diagrams with a tadpole-type 
self-energy insertion (the loop consisting of a single closed photon line) since it is set to zero 
in dimensional regularization. Thus we are left with 24 one-photon loop diagrams which we 
group into the (separately crossing-symmetric) classes I-IX, shown in Figs. 2-10. Due to an 
increasing number of internal pion propagators their evaluation rises in complexity. In order 
to simplify all calculations, we employ the Feynman gauge, where the photon propagator is 
directly proportional to the Minkowski metric tensor g'^'^. We can now enumerate the analytical 
expressions for A{s, u) and B{s, u) as they emerge from the 9 classes of contributing one-photon 
loop diagrams. 

Class I : For this class of diagrams the tree level amplitude in eq.(4) gets multiplied with the 
wavefunction renormalization factor of the pion: 



2a 



7r{u — 1) 



(6) 



Class II : These diagrams involve the once-subtracted (off-shell) self-energy of the pion: 



B{s,u) 



(11) 



a 



^{u — 1) 



"U + 1 

2^1/1/ + 2 — In(l-M) 



u 



Class III : These diagrams generate a constant vertex correction factor: 



a 



tt{u 



36 



uv 



Class IV : These diagrams generate a w-dependent vertex correction factor: 



i?(s,w)(^v) 



a 



tt{u — 1) 



1 u — 1 



(7) 



(8) 



(9) 
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Class: III 




Figure 4: One-photon loop diagrams of class III. 



Class: IV 




Figure 5: One-photon loop diagrams of class IV. 



Class: V 



Figure 6: One-photon loop diagrams of class V. 



Class: VI 



to, 



Figure 7: One-photon loop diagrams of class VI. 
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Class: VII 



Figure 8: One-photon loop diagrams of class VII. 



Class: VIII 




Figure 9: One-photon loop diagrams of class VIII. 



Class V : These diagrams generate also a ^-dependent vertex correction factor: 



B{s,u) 



(V) 



a 



47r{u - 1^ 



- 4^uv + 5 + - + ln(l - u) 



u 



Class VI : These diagrams generate a u- and Q-dependent vertex correction factor: 



B(s, u 



,{vi) 



—, -T -4et7V' + 5 + - + ^ln(l-M) 

-liu + 3 + 2g) Cdx f\ 

Jo Jo K[u,Q) 



(10) 



(11) 



with the cubic polynomial K{u,Q) = y + {1 — u)x{l — y) + Qx{l — x)y in two Feynman 
parameters x^y. Note that in the physical region m < 1,Q > the denominator K{u,Q) > 
is strictly positive and therefore the integral in eq.(ll) is well-defined and can be performed 
easily numerically. Interestingly, for the classes I- VI considered up to here, the contributions 
to the amplitude A{s,u) vanish identically in each case. 
Class VII : These irreducible diagrams give: 



a 



5 — Sm , , , 1 
ln(l - m) - 2 



An(u — 1) L u 
a ( 1 + u 1 — 3u 



u - 1 



2Li2(M) - - 



u 



ln(l - m) + 



1 



u-\ 



vr 



- - 2Li2(^i) 



(12) 



(13) 



47r(M - 1) 

where \j\-2iu) = Y^'^=i^~'^W^ = u [x — u)] Inx denotes the conventional dilogarithmic 

function. 

Class VIII : These irreducible diagrams give: 



(14) 



Class: IX 





Figure 10: One-photon loop diagrams of class IX. 



B(s, u 



,(viii) 



- Cdx Cdv ^^^^ ~ 
27r7o io K{u,Q) 



(15) 



The (singular) integral in eq.(14) involving l/K{s,Q) requires a precise prescription for its 
(numerical) treatment. First, s > 1 is given an infinitesimally positive imaginary part s + zO"*", 
which makes the integral complex-valued. In order to compute the only relevant real part of 
it, one performs analytically the Feynman parameter integral /J^ dy and converts the occurring 
logarithm ln[a;(l — s)] into ln[a;(s — 1)]. The remaining integral dx is free of poles and therefore 
unproblematic for a numerical treatment. 

Class IX : These box diagrams are most tedious to evaluate and they give: 



A{s, u 



I (IX) ^ ^ 



1 



27r - 1 
1 



u + l 
2u{u - 1) 



ln(l -m) 



1 



y -Li2(M) 



+ 



2QJA + QL{Q) 



Q + t-At L\-t) - 4Q L'^iQ) + 2Q^A^t L{-t) 
+ I dx dy 



1 r(l_a;)y(l-y) 











Kis,Q) 



{2t - A)xy 



/ [1 — x){l + xy — y — xyz) x{l — xy){l — z) 
/o \ G'2(s,t,Q) ^ G'2(M,t,Q) 



where we have introduced the frequently occurring logarithmic loop function: 



1 V4 + X + 
L(x) = ^ln 



(16) 



(17) 



and the quartic polynomial G{s^ t,Q) = y + [1 — s)x{l — y)z — tx{l — x)y{l — z) + Qx{l — x)yz 
in three Feynman parameters x, y, z. The other contribution reads: 



a 



27i{u - 1) 
2t -4 



+ 



2u^ + U-1 
2^^ 

4(e/i? + ln(l 



ln(l -m) - 



u + l 
2u 



+ 



IT 



M — 1 V 6 



tL(-t) +Li2(w) -Li2(l 



w] 



In w 



--ln\l-w) + U2{h 



Li2(/^^ 



Vt^ - 4t / dx 
Jo 



\n[l + Qx{l - x)] 



u - 1 



2QJA + QL{Q) 



1 — tx{l — x) 
Q-t + At L\-t) + AQ L\Q) - 2QVA^ L{-t) 
1 ^(l-xy){y-l) 



+ (u 



I dx dy 
Jo 



K{u,Q) 



_ . f x{l - x)y{z - 1) 1 - xy + x^y + xz - x'^yz 

-[zt ^)xyj^azy G%s,t,Q) + G^{u,t,Q) 



[18) 



with the auxihary variables: 



^=^(^1-^4^), h^ = l(t±VW^), (19) 

depending ont = 2 — s — u — Q. In the physical region u < l,t < 0,Q > the denominator 
polynomial G{u, t,Q) > is strictly positive and therefore the integrals involving t, Q) 

are well-defined and suited for numerical evaluation. The (highly singular) integrals involving 
1/G^(s, t, Q) with s > 1 (after the substitution s ^ s + i 0^ which shifts the double-pole off the 
real integration contour) must however be treated partly analytically. The (solvable) integrals 
dz over the third Feynman parameter z appear in the following two versions: 

lo '^^ G^{s,t,Q) " [l-tx{l-x)]K{s,Q) ' 

\ x^z X 1 , K(s,Q) 
(j^ — I In \ /oil 

G'is.t.Q) N{s,Q + t)K{s,Q) N^{s,Q + t) [1 - tx{l - x)]y ' ^ ' 

giving rise to a new denominator polynomial N{s, Q + t) — {s — — y) — {Q + t){l — x)y > 
which is strictly positive for s > 1 and t < —Q. Although the latter condition t < — Q 
does not hold in the entire physical region, it represents (fortunately) no restriction for the 
(experimentally selected) events in the Coulomb peak, where Q < 0.1 (or even less). By 
analyzing the condition t < —Q in the laboratory frame one deduces the inequality for the 
ratio of photon to pion beam energy, uj/Ei> [^Q{^ -\- Q) — Q\l2 ~ \fQ — Ql2., which is indeed 
well satisfied for the experimentally selected events. In the next step, the terms in eqs. (20,21) 
with a factor 1/A'(s, Q) are integrated once more analytically, /g^ rfy, and for isolating the only 
relevant real parts the emerging logarithms are converted into logarithms of absolute values. 
The same conversion into real- valued logarithms is done with the term lni^(s, Q) in eq.(21), 
this time without solving the integral /q^ dy analytically. The remaining single and double 
integrals are free of poles and thus suitable for a straightforward numerical evaluation. The 
same strategy for evaluating the occurring Feynman parameter integrals applies to the crossed 
amplitudes, A(u^s) and B{u,s), with the variables s and u interchanged. An alternative to 
the Feynman parameter representation of complicated loop functions in eqs. (11,14-18) is the 
(dispersive) spectral function representation based on the Kramers-Kronig dispersion relation 
(see section 3 in ref.[ll]). However, at non-zero Q the corresponding spectral functions (or 
imaginary parts) entering these principal value integrals turn out to be very cumbersome. If an 
expansion up to linear order in Q suffices, one can do further steps in the analytical integration. 

As a good check one verifies that the ultraviolet divergent terms proportional to ^uv cancel 
in the total sum for B{s,u). Interestingly, the contributions to the other amplitude A{s,u) 
written in eqs. (12, 14, 16) are ultraviolet and infrared finite in each case. 
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4 Infrared finiteness 

In the next step we have to consider the infrared divergent terms proportional to ^/i?. Inspection 
of eqs. (6,18) reveals that they scale with the (tree level) Born amplitude written in eq.(4). As 
a consequence of that feature, the infrared divergent virtual (loop) corrections multiply the 
point-like differential cross section da^^^^ /dcodQ^dQ-,,^ by a (s <-> crossing-symmetric) factor: 

Cm. (22) 



virt 



2a 

TT 



2f — 4 
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The (unphysical) infrared divergence ^jff gets canceled at the level of the (measurable) cross 
section by the contributions of soft photon bremsstrahlung. In its final effect, the (single) soft 
photon radiation multiplies the tree level cross section da^'^^y dLodQ^dflT^ by a (universal) factor 
[11, 15]: 



4oft = a 



4-d 



\\<\ 



(23) 



which depends on a small energy cut-off A. Working out this momentum space integral by the 
method of dimensional regularization (with d > 4) one finds that the infrared divergent correc- 
tion factor 5^;^^ ~ ^7ij in eq.(22) gets eliminated and the following finite radiative correction 
factor remains: 



r(cm) 
real 



2 + 



At 



s + l + Q 



L{-t) 



s + 1 

In — - H ; 

2A 2 s 



1 



In s 



^(s - 1 + QY + 4g 



In 



S + 1 + Q + J{s-l + Qf + AQ 



+(i-2) l^dx 
Jo 



1 + Qx 



In- 



2^ 

1 + Qx 



W 



2W[l-tx{l-x)] s + l + Qx-W 



(24) 



with the abbreviation W — ^'(s -|- 1 -|- Qx)'^ — 4s -|- Astx{l — x). We note that the terms beyond 
those proportional to ln(m,r/2A) are specific for the evaluation of the soft photon correction 
factor (5soft in the 7r~7* center-of-mass frame with A an infrared cut-off therein. As an aside we 
mention that at Q = and in backward direction, t = —{s — l)'/s, the integral in eq.(24) can 
be solved and one obtains the following handy formula: 



r(back) 
real 



a 



TT S2 _ 1 



l-{s^ + l)]ns In^ 



2A 



■ln's + 2Li2(l -s) 



(25) 



We can now present some numerical results for the radiative corrections to the pion-nucleus 

bremsstrahlung reaction 7r~Z — >■ 7r~Z'y. The radiative correction factor 6^^^^^ due to soft photon 
emission vanishes for t = Q = at all s > 1. The full curves in Fig. 11 show its values for 
three selected n^j center-of-mass energies \/sm-„ = (2, 3, 4)m,r as a function of Q (the photon 
virtuality divided by m^) choosing the extremal momentum transfer tmim where — 2stinin = 

\2 



(s — 1) + Q{s + 1) + (s — l)y (s -f 1 + Qy — As corresponds to the backward direction in the 
7r~7* center-of-mass frame. The soft photon detection threshold has been set to A = 5 MeV. 
One observes negative values up to —3% with a weak dependence on Q (in particular in the 
region of the Coulomb peak Q < 0.1). The three dashed lines in Fig. 11 show for comparison the 
case of perpendicular scattering with momentum transfer t± reduced to about one half of tmin, 
where —2st± = (s — l)' + (5(s + l). At higher center-of-mass energies the excitation of the broad 
p(770)-resonance becomes prominent, = mp/m.jr — 5.5. However, this kinematical region is 
not of interest for the extraction of the pion polarizabilities. The numbers in Fig. 11 should be 
considered only as indicative since there are yet the (finite) virtual radiative corrections due to 
the photon-loops. A comparison with the case of real pion Compton scattering [11] suggests that 
these will reduce the magnitude of the radiative corrections to pion-nucleus bremsstrahlung. 
At the present stage, a proper quantification of these (finite) virtual radiative corrections is not 
possible since that requires specification of the actual experimental conditions, such as pion 
beam energy, constraints on detectable energy and angular ranges, etc. However, their full 
implementation into the analysis of the COMPASS data is currently planned. 
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Figure 11: Radiative correction factor S^.^^^ for virtual pion Compton scattering 7r~7* — > 7r~7 
as a function of the photon virtuahty Qm^. The soft photon detection threshold has been set 
to A = 5 MeV in the n^j center-of-mass frame. Full and dashed lines correspond to backward 
and perpendicular scattering, respectively. 



5 Inclusion of pion structure via polarizability difference 

In the same way as Akhundov et al. [13], we have treated so far the pion as a structureless spin- 
boson in our calculation of the radiative corrections. Now we go beyond this approximation 
which is valid only near the 7r^7 threshold. The leading pion structure relevant for (virtual) 
Compton scattering is given by the difference of its electric and magnetic polarizability — /J^r- 
In an effective field theory approach it is easily accounted for by a new two-photon contact- 
vertex proportional to the squared electromagnetic field strength tensor, F^^F^^ — 2(B^ — E^). 
The S-matrix insertion following from this higher-order (gauge-invariant) effective 7r7r77- vertex 
reads: 

87ri/3^m^(^A;i ■ k2 ei ■ — ■ k2 ■ kij , (26) 

where one photon {ki,ei) is ingoing and the other one (A:2,e2) outgoing. At tree level (see 
left diagram in Fig. 12) the polarizability vertex eq.(26) gives in the kinematical situation of 
7r~7o — > 7r~7 rise to the (constant) contributions: 

A{s, u) (P°i^) ^-B{s,u) (P°'^) = - , (27) 

a 

to the amplitudes A(s,u) and B{s,u) parameterizing the laboratory T-matrix Tiab- In order 
to prevent any misunderstandings we stress that when writing (merely) /^TrmTr for the coupling 
strength in eq.(26), an electric and a magnetic pion polarizability, equal in magnitude and 
opposite in sign = —^t^, are always both included. 

We reinterpret now the one-photon loop diagrams of section 3 in such a way that the two- 
photon contact- vertex represents one polarizability vertex proportional to /^Trm^. In order to 
distinguish it from the leading order contact- vertex Snia ei • 63 of scalar QED we have symbol- 
ized it in Figs. 12-15 by an open square. Going through the classes F- V and evaluating the 
loop diagrams with the S-matrix insertion from the polarizability vertex we find the following 
contributions to the amplitudes A{s,u) and B{s,u). 
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Class: r 



Figure 12: The open square symbolizes the electromagnetic interaction proportional to the 
(dominant) pion polarizability difference aTr—lS-jr. Shown are the tree diagram and the diagrams 
of class r involving the pion wavefunction renormalization factor. 



Class: ir 




Figure 13: One-photon loop diagrams of class IF including pion structure via its polarizability 
difference. 
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Figure 14: One-photon loop diagrams of class IIF 
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Class: IV' 



V 



Figure 15: One-photon loop diagrams of classes IV' and V. 
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\nK{s,Q) 



(32) 



with the polynomial K{u,Q) defined below eq.(ll). The strategy to compute the integrals 
involving l/K{s,Q) and \aK{s,Q) has been described in section 3. 
Class IV': 



,{IV') 
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(33) 
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U2{u) - — 
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X 



46i?v^ L{-t) + Li2(u;) - Li2(l - w) + w - w) 



(34) 



(35) 



with the auxiliary variable w defined in eq.(19). The contributions of the vertex type diagrams 
(classes III and IV in Figs. 4,5) vanish due to some special features of the polarizability vertex 
eq.(26). 

At first count the ultraviolet divergent terms proportional to ^jjv do not drop out in the 
total sums for A{s,u) and B{s,u). They disappear however after interpreting the coupling 
constant in eq.(27) as a bare one and splitting it as P^^"^^^ = /?7r(l — « ^uv/2t:) into the physical 
coupling constant (3^^ and a counterterm piece. Note that the same renormalization procedure 
has rendered the radiative corrections to real pion Compton scattering ultraviolet finite [11]. 
The infrared divergent terms proportional to E,ir showing up in eqs. (28,35) are again canceled 
by the soft photon bremsstrahlung contributions. As a result one has to multiply the partial 
differential cross section da'^^°^'' / dijjdVL.ydVL.^ linear in which arises from the interference of the 
Born and polarizability terms in eqs. (4,27), with the radiative correction factor 6^^™^ written 
down in eq.(24). In the case of real pion Compton scattering it has been observed in ref. [11] 
that the inclusion the leading pion structure effect (via — ~ —2/?^) does not change the 
relative size and angular dependence of the radiative corrections. Due to the dominance of real 
over the virtual radiative corrections (for small enough infrared cut-off A) one can expect that 
the same feature holds for the pion-nucleus bremsstrahlung reaction tt~ Z ti~ Z'^. 
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z 




Figure 16: Two-photon exchange diagrams specific for the contact-vertex of scalar QED. The 
double-hne symbohzes the heavy nucleus of charge Z. 

Appendix: Specific two-photon exchange contributions 

For high nuclear charges Z the one-photon exchange approximation is insufficient for analyzing 
experimental data of pion- nucleus bremsstrahlung tt~Z —>■ n^Z'-f. A non-perturbative treat- 
ment of the scattering of the charged particles in the strong nuclear Coulomb field is required 
in this case. Integrating the tree-level cross section da^^^^ / dujctVL^dVL.,^ of scalar QED over the 
angles (6'i, 6*2,0), and adapting the result of refs.[16, 17] obtained for electrons to pions, one 
may write the photon spectrum at high beam energies Ei in the form: 



with uo = El — E2 the photon laboratory energy. The function ^!^{Za) represents the non- 
perturbative corrections beyond the one-photon exchange. For electrons this function has been 
calculated as '^e{Za) = — Z"^ Yl'^=i[nij^^ + Z'^a^)]~^ in refs.[16, 17]. In the case of a lead 
or a nickel target {Z = 82 or 28) the negative infinite series "^eiZa) amounts to —0.3316 or 
—0.0484, i.e. about a —5% or —0.7% correction compared to the leading logarithm plus constant 
in eq.(36). At this point it would be most important to clarify, whether these non-perturbative 
corrections are the same for pions and electrons or not. 

As a supplement to the subject, we evaluate in this appendix those two-photon exchange 
contributions (of order Z^ in the transition amplitude) which are specific to the two-photon 
contact- vertex of scalar QED. The pertinent one-loop diagrams are shown in Fig. 16. Addition- 
ally possible box and crossed box type diagrams, where one exchanged virtual photon and the 
emitted real photon are connected by the contact- vertex vanish as a consequence of the gauge 
condition eo = 0. For the evaluation of the (triangular) loop we use a non-relativistic (static) 
propagator i/{lQ+i O"*") for the intermediate heavy nucleus. By dividing out the factor — eZ/|g p 
we can translate the two-photon exchange processes into a correction to the T-matrix Tiab of 
virtual pion Compton scattering. Putting all pieces together we find from the photon-loop 
diagrams in Fig. 16 the following contribution: 



It has a structure very similar to the T-matrix in Born approximation (see eqs.(3,4)). Making 
this comparison explicitly, one deduces that the coupling terms e* ■ pi^2 in eq.(37) carry the 
relative suppression factors ^laZm^,^/Q /4:E2,i, respectively. Although naZ/A ~ 0.47 is already 
sizeable for a lead target {Z = 82), the ratio \q\/Ei^2 of (virtual photon) momentum transfer 
to beam or scattered pion energy reduces it by at least three orders of magnitude. We are 
considering here the typical kinematical situation, £^12 > 50 GeV and \q\ < 50MeV, where 




(36) 




(37) 
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high-energy pion-nucleus bremsstrahlung vr^Z ^ is used to extract the pion polarizabil- 

ities. In passing we note that the diagrams in Fig. 16 with the contact-vertex replaced by the 
polarizabihty vertex give rise to a T- matrix that differs from eq.(37) merely by an additional 
prefactor P-^ml.Q/2a ~ — 7 - 10~^ Q (using = — ~ 3- lO^^fm^). We can therefore conclude 
that the two-photon exchange mechanisms in Fig. 16, specific to the pseudo-scalar pion, are 
negligibly small. 
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